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FOREWORD 


This  final  technical  report  was  prepared  by  the  Massachusetts  Institute  of 
Technology,  Cambridge,  Massachusetts,  on  Contract  AF  33(616) -6076  for  the 
Aeronautical  Research  Laboratory,  Office  of  Aerospace  Research,  United  States 
Air  Force.  The  work  reported  herein  was  accomplished  on  Task  70138,  "Research 
on  a  Novel  Technique  of  Measuring  Very  High  Temperature  Gas"  of  Project  7064, 
"Research  on  Aerodynamic  Fields",  under  the  technical  cognizance  of  Dr.  Max 
G.  Scherberg  and  Mr.  Erich  Soehngen  of  the  Thermo -Mechanics  Research  Branch 
of  ARL. 

Part  I  is  based  on  a  Mechanical  Engineering  thesis  entitled,  "An  Analytical 
Study  of  the  Effect  of  Vibrations  on  Heat  Transfer  from  a  Heated  Horizontal 
Cylinder"  presented  in  I960  by  T.  Chiang,  Research  AssistcUit,  Research  Lab¬ 
oratory  for  Heat  Transfer  in  Electronics,  MIT. 

Appreciation  is  extended  to  Prof.  R.  J.  Nickerson  for  his  valuable  original 
suggestions  and  for  his  critical  reading  of  the  final  manuscript.  Prof.  Nickerson 
suggested  the  analytical  approach  by  means  of  which  the  free-convection  solution 
described  in  Mr.  Chiang's  thesis  was  obtained. 
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ABSTRACT 


In  this  report,  the  basic  equations  and  the  boundary  condi¬ 
tions  which  govern  the  problem  of  coupled  transverse  vibrations  and 
free  convection  from  a  heated  horizontal  cylinder  are  presented.  By 
applying  a  method  developed  by  C.  C.  Lin,  it  is  shown:  1)  that  the 
presence  of  harmonic  oscillations  modify  the  steady-flow  solution 
only  when  pressure  gradients  are  present;  2)  that  the  modifying  forces 
have  their  most  pronoianced  effect  on  the  fluid  closest  to  the  surface; 
and  3)  that  the  product  a<i>  is  a  measure  of  the  magnitude  of  the  modify¬ 
ing  forces.  The  use  of  the  quantity  emj  as  a  measure  of  the  magnitude 
of  the  influence  of  vibrations  on  free  convection  agrees  with  experi¬ 
mental  correlations. 

By  transforming  the  differential  equations  into  dimensionless 
form,  it  is  shown  that  four  dimensionless  parameters  are  needed  to  fully 
describe  the  flow.  Four  sets  of  dimensionless  parameters  and  their 
corresponding  differential  equations  are  presented;  the  advantages  and 
disadvantages  of  each  set  are  discussed.  A  perturbation  method  is  applied 
to  one  set  of  equations,  and  the  zeroth-order  solution  is  obtained.  This 
zeroth-order  solution,  which  corresponds  to  free  convection,  agrees  with 
Hermann's  analysis  for  a  heated  horizontal  cylinder. 
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I.  INTRODUCTION 


During  the  past  few  years >  considerable  Interest  has  arisen 
with  regard  to  the  Interaction  between  vibrations  and  convective  heat 
transfer.  These  coupled  phenomena  have  been  studied  experimentally  by 
two  methods.  The  first  method  utilizes  a  sound  field  In  a  fluid  medium 
surrounding  a  heated  object;  In  the  second  method,  a  heated  object  Is 
caused  to  vibrate  by  mechanical  means  In  an  otherwise  undisturbed  fluid 
medium.  Both  methods  result  In  a  relative  velocity  vector  between  the 
heated  surface  and  the  fluid  medium  which  changes  periodically  In  magni¬ 
tude  6Lnd  direction. 

One  of  the  earliest  Investigations  of  the  Influence  of  vibra¬ 
tions  on  convective  heat  tremsfer  was  performed  by  Martlnelll  and 
Boelter^^^  using  an  apparatus  In  which  a  heated  horizontal  cylinder  was 
vibrated  vertically  In  a  tank  of  water.  For  sxifflclently  Intense  vibra¬ 
tions,  Martlnelll  and  Boelter  reported  that  the  overall  heat- transfer 

coefficient  was  Increased  by  as  much  as  four  times  Its  value  without 

(2  3) 

vibrations.  More  recently,  Feuid  and  Kaye'  *  have  made  experimental 
Investigations  of  the  Influence  of  acoustical  and  mechanical  vibrations 
on  heat  tremsfer  by  free  convection  from  horizontal  cylinders.  Their 
data  shov  that  when  the  Intensity  of  vibrations,  defined  as  the  product 
of  amplitude  and  frequency,  af,  exceeds  a  certain  "critical  value" 
(approximately  0.3  ft/sec),  the  coefficient  of  heat  transfer  Increases 
significantly.  The  Influence  of  vibrations  on  heat  tremsfer  has  also 
been  studied  by  Lemllch^^^,  Tsul^^\  Shlne^'^\  and  Sprott,  Holman,  and 
Durand . 
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(3  a) 

Fand  and  Kaye'^*  have  carried  out  flow-visualization  studies 
which  demonstrate  that  the  boundary  layer  flow  about  a  heated  horizontal 
cylinder  changes  radically  In  the  presence  of  super-critical  vibrations. 
Their  photographs  of  boundary-layer  flow,  using  smoke  as  the  Indicating 
medium,  show  that  vertical  (transverse)  vibrations  affect  the  boundary  layer 
in  a  different  way  than  do  horizontal  (transverse)  vibrations.  In  the 
case  of  horizontal  super-critical  vibrations,  a  vortex-type  flow,  called 
thermoacoustic  streaming,  was  observed  to  occur  about  the  cylinder;  whereas, 
in  the  case  of  vertical  super-critical  vibrations,  a  wide  zone  of  turbu¬ 
lence  surro\inded  the  test  cylinder.  Thermoacoustic  streaming  has  also 

(a) 

beun  observed  by  Sprott,  HoLnan,  and  Durand'  '  by  shadowgraph  methods. 

The  analytical  aspects  of  the  coupling  between  vibrations  and 
free  convection  constitute  a  very  difficult  problem — to  date  no  theoreti¬ 
cal  analysis  of  this  problem  has  appeared.  However,  analyses  of  several 

(Q\ 

related  problems  have  been  successfully  carried  out.  Stokes  ,  and 
later  Rayleigh^ analyzed  the  case  of  an  infinitely  long  plate  execut¬ 
ing  harmonic  oscillations  parallel  to  its  own  plane  In  a  viscous  fluid 
otherwise  at  rest.  This  case  was  relatively  simple  because  there  was  no 
temperature  difference  between  the  plate  and  the  surrounding  fluid,  and 
also  because  the  motion,  being  Independent  of  the  coordinate  along  the 
plate,  satisfies  a  linear  differential  equation.  Stokes'  solution  is 
known  as  "shear-wave"  flow.  Schlichting^^^^  considered  analytically 
the  case  of  a  cylinder  executing  harmonic  oscillations  In  a  fluid  at 
rest  with  no  temperature  difference  between  the  cylinder  and  the  fluid. 

He  used  boundary-layer  approximations  and  then  applied  an  Iterative 
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approach  against  a  steady  mean  flow,  which  Is  Identically  sero.  In  order 
to  solve  the  governing  eqxiatlons.  He  found  that  the  first  approxima¬ 
tion  to  the  solution  represents  a  periodic  flow,  and  the  second  approxima¬ 
tion  to  the  solution  contains  a  steady-flow  tern  which  Is  Independent  of 
time.  Schllchtlng's  analytical  results  explained  the  acoustic  stream- 
Ing  phenomena  observed  earlier  by  Andrade'  .  Schllchtlng's  original 
analysis  of  acoustic  streaming  has  bean  refined  and  extended  by  Holtsmark 
et  al^^^\  Raney  et  al^^\  and  Skavlem  and  Tjotta^^^^.  The  effect  of 
free-stream  oscillations  on  laminar  boundary  layers  have  been  considered 
by  Ostrach^^^\  Moore^^'^\  emd  Cheng  and  Elliott^^^^ . 

(19) 

Llghthill'  '  studied  the  laminar  boundary-layer  response  to 
fluctuations  in  stream  velocity  by  using  an  approximate  boundary-layer 
method,  Nlckerson^^*^^ ,  in  a  combined  theoretical  and  experimental 
Investigation,  considered  the  effect  of  free-stream  oscillations  on  the 
laminar  boundary  layer  from  a  flat  plate  amd  the  associated  forced-con¬ 
vective  heat-transfer  problem;  the  fluctuations  In  the  free  stream  were 
assumed  to  be  small  compared  to  the  steady  mean  velocity.  Nickerson's 
analysis  predicted  a  negligible  change  in  the  coefficient  of  heat-trans¬ 
fer  from  the  plate  surface;  experimental  data  corraborated  this  analyti-* 
(21) 

cal  result.  Lin  made  a  theoretical  analysis  on  the  motion  in  boundary 
layers  with  a  rapidly  oscillating  external  flow.  His  analysis  was  not 
limited  to  small  oscillations  In  the  free  stream.  Lin's  results  Indi¬ 
cated  that  the  Influence  of  oscillatory  motion  on  the  mean  boundary- 
layer  profile  Is  contingent  upon  the  existence  of  pressure  gradients 
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outside  the  boimdary  layer  and  predicted  the  possibility  of  a  periodic 
back  flow. 

The  foregoing  brief  literature  review  shows  that  several 
experimental  Investigations  of  the  interaction  between  vibrations  and 
free  convection  have  been  performed,  but  no  theoretical  analysis  of  this 
Interaction  has  bean  carried  out.  This  report  represents  a  first  step 
toward  the  solution  of  this  complex  problem.  More  specifically,  this 
report  contains  a  formulation,  partial  solution,  and  discussion  of  the 
system  of  differential  equations  which  describe  the  coupling  between 
free  convection  and  transverse  vibrations  of  a  heated  horizontal  cylinder. 
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II.  BASIC  EQUATIONS  AND  THEIR  PROPERTIES 

2.1  Coordinate  System 

The  system  of  differential  equations  for  the  velocity  and 
temperature  distributions  near  a  heated  horizontal  cylinder  subjected 
to  transverse  harmonic  vibrations  will  be  derived  in  this  report.  It 
will  be  assumed  that  the  cylinder  is  sufficiently  long  so  that  there 
is  no  flow  in  the  axial  direction;  hence  the  problem  is  two-dimensional. 

An  oscillating  reference  frame  which  is  fixed  to  the  cylinder 
will  be  used  for  the  analysis;  thus^  for  an  observer  standing  on  the 
cylinder,  the  cylinder  appears  stationary,  and  the  fluid  at  infinity 
executes  harmonic  oscillations.  This  oscillating  frame  of  reference 
will  give  results  which  are  equivalent  to  those  using  a  reference  frame 
fixed  in  space,  due  to  the  fact  that  any  inertia  force  per  unit  volume 
is  balanced  by  an  equal  pressure  gradient.  This  point  has  been  explained 
by  Llghthill^^^^ 

A  curvilinear  orthogonal  system  of  coordinates  will  be  used 
to  simplify  the  analysis:  The  x-axis  of  this  coordinate  system  Is  along 
the  wall  of  the  cylinder,  and  the  y-axis  is  perpendicular  to  it,  as  shown 
in  Fig.  1. 

2.2  Basic  Equations 

(22) 

Tollmlen  derived  expressions  for  the  complete  Navler-Stokee 
aquations  in  the  coordinate  system  indicated  in  Fig.  1.  These  equations 
will  be  used  in  a  modified  form  in  the  analysis  which  follows.  In  this 
modification,  the  density  variation  will  be  Introduced  exclusively  in 
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the  buoyancy  term  in  the  momentm  equatlone*  and  the  additional  aseump- 
tion  will  be  made  that  the  radius  of  the  cylinder  is  l€u*ge  compared  to 
the  boiindary-layer  thickness. 

By  applying  the  usual  boundary-layer  assumptions  and  adding 
the  buoyancy-force  term  in  the  momentum  equations  derived  by  Tollmlen, 
the  bo\mdary-layer  equations  become: 


Continuity: 


^  +  jLj:  = 

dx 


0 


Momentiun: 


Energy: 


9  t 


+  u 


S  X 


+  V 


3G 

S>y 


(1) 

(2) 

(3) 


It  should  be  noted  that  these  equations  do  not  apply  near  the  upper  por¬ 
tion  of  the  cylinder  where  the  normal  component  of  velocity  becomes  of 
the  same  order  of  magnitude  as  the  tangential  component,  because  this 
is  contrary  to  the  boundajry-layer  assumptions. 

Outside  the  boundary  layer,  the  nonsteady  Bernoulli  equation 

holds. 


(>  3x  3t  “3x 


(4) 


2.3  Boundary  Conditions 

The  most  general  form  of  the  boundary  conditions  for  Eqs.  (l) 
through  (4)  are: 

when  y  =  0,  u  =  v  =  0,  0=  0 


when  y-^eo,  u-*U(x,  t) ,  O  ^0 


(5) 
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K- 


ilj 

t 


These  boundary  conditions  can  be  made  more  specific  by  choosing  vari¬ 
ous  types  of  velocity  conditions  far  from  the  cylinder.  The  class  of 
boxuidary  conditions  considered  in  this  report  are  those  in  which  the 
fluid  far  from  the  cylinder  executes  harmonic  vibrations  in  one  direc¬ 
tion.  The  mathematical  expression  for  this  class  of  boiindary  condi¬ 
tions  is: 

when  y-#  00)  u-*U  =  2  acj  sin  Ot  sin  +  X)  (6) 

K 


This  boundaiy  condition  corresponds  to  the  cylinder  being 
vibrated  in  an  otherwise  undisturbed  fluid  with  the  coordinate  system 
fixed  on  the  cylinder.  It  also  approximates  the  case  of  a  stationary 
cylinder  in  a  strong  transverse  so\md  field  when  the  wave  length  of 
the  sound  field  is  large  compared  to  the  diameter  of  the  cylinder.  The 
direction  of  these  vibrations  can  be  specified  by  an  appropriate  choice 
of  A  in  Eq.  (6).  Thus,  for  vertical  vibrations,  X=  Oj  for  horizontal 
vibrations,  X=  it/2. 

2.4  Method  of  Solution  Developed  by  C.  C.  Lin 

(21) 

C.  C.  Lin  has  developed  a  method  of  attacking  unsteady 
forced-flow  boundary-layer  problems  which  supplies  a  great  amount  of 
physical  insight.  This  method  applies  to  periodic  flows  when  the  ’’AC" 
boundary-layer  thickness  is  small  compared  to  the  steady  boundary-layer 
thickness;  that  is,  when 


0 

T~f 


« 1 


(7) 
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In  Lin's  method,  the  assumption  is  made  that  the  flow  can 
be  divided  into  a  steady  part  and  an  oscillating  part  in  the  follow¬ 


ing  form; 

u(x»  y»  t)  =  u(x,  y)  +  Ui(x,  y,  t)  Uj  =  0 

v{x,  y,  t)  =  v(x,  y)  +  Vi(x,  y,  t)  7^  =  0 

U(x,  t)  =  U(x)  +  Ui(x,  t)  Ui  =  0 

p(x,  t)  =  p(x)  +  pi(x,  t)  Pi  =  0 

d  (x,  y,  t)  =  Q(x,  y)  +  di(x,  y,  t)  Oj  =  0 


(8) 


With  these  assumptions,  the  boundary-layer  equations  are 
separated  into  two  sets:  one  fcr  the  steady  part  of  the  flow  and  one 
for  the  oscillating  part.  The  steady  set  of  equations  is  obtained  by 
taking  the  time  average  of  the  complete  equations.  The  equations  for 
the  oscillating  flow  are  obtained  by  subtracting  the  steady  equations 
from  the  complete  equations  according  to  the  following  procedure. 

First,  Eq.  (4)  is  used  to  eliminate  the  pressure  term  from 
Eq.  (2).  Then  Eqs.  (8)  are  substituted  into  Eqs.  (l),  (2),  and  (3). 
This  manipulation  results  in  the  following  set  of  equations: 


ui)  +  ^i)  =  0 

5  y 


(9) 


3(u  +  Ui)  _  0(u  +  Uj)  5  Ui 


^  (U  +  Ui) 


+  (u  +  Ui)  - -  +  (v  +  Vj) 

^  (U  +  Ui)  .  ^^(u  +  Ui) 


5  y  ~  3  t 


(10) 


d  X 


^  Vu  •  uj/ 

V  ^  y2  +  +  ®i)  Sin  I 
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(11) 


dSi  0(6+  6i)  0(6  +  di) 

JL- 

f®p  •)  y* 


< 


When  the  time  averages  of  Eqs.  (9),  (lO),  and  (ll)  are  taken,  the  follow¬ 
ing  expressions  for  the  time  independent  part  of  the  flow  are  obtained: 


0  u 
Ox 


(12) 


•)u,-iu_A 

-=—  +  V  -5^  -  V 

Oy 


^5  +  g/?  6 


sin  I  + 


(13) 


-  36  ^ -  06  _  k  326  ^ 


(14) 


where 


(15) 


2  06i  Odi 

Ox  ^  0  y 


(16) 


B.C.  y=0,  u  =  v  =  0,  6  = 

y+*0,  u 0,  6  "+0 


(17) 


The  equations  describing  the  oscillating  part  of  the  flow  are 
obtained  by  subtracting  Eqs.  (12),  (13),  and  (14)  from  Eqs.  (9),  (10), 
and  (11).  These  equations  are: 


Ouj  Ovi  ^ 

O  x  O  y  ~  ° 


(18) 
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(19) 


■  ~s~P 


5Ui 

at 


-Jl  -  fc 

at  " ^ Cp  Tr  ~  " 


>  .  aui  jui 

where  /  =  -^  -  u^  -  v,  -/+ 

_  ,  d'li 

-u  - - -  -r —  + 

dx  ay 

-“‘  TT  -  ’i  51 


> 


+ 


.  dOi  . 

-  u  ^'-  -  V  -r -  + 

a  X  Oy 

dO  oo 

-  aT  ‘  oT 


B.C. 


y  =  0,  Ui  =  ▼!  =  0, 


=  0 


y-^e«  ,  Ui-^2a<i>  sin  wt  sin  (^  +  X) ,  ^  0 


(20) 


(21) 


(22) 


(23) 


Lin  suggests  an  iterative  approach  when  his  method  is  applied 
to  forced  flows.  In  forced  flows >  the  coupling  between  the  velocity  and 
temperature  fields  can  be  considered  negligible  in  many  cases.  Hence, 
for  forced  flows,  the  equations  can  be  solved  for  the  velocity  distribu¬ 
tion  independently  of  the  temperature  distribution.  However,  the  strong 
coupling  between  the  velocity  and  the  t«aperature  fields  which  exists  in 
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the  case  of  free  convection  presents  a  serious  obstacle  to  obtaining 
a  solution  by  this  approach;  for,  with  free  convection,  the  velocity 
and  temperature  distributions  cannot  be  solved  independently. 

2.5  Characteristics  of  the  Flow  Shown  by  Lin's  Method 

Although  Eqs.  (12)  through  (23)  are  difficult  to  solve, 
they  do  show  several  interesting  characteristics  of  luisteady  flows. 

First,  they  show  that  the  steady  part  of  the  flow,  when  oscillations 
are  present,  is  different  from  the  steady  solution  with  no  oscillations— 
due  to  the  terms  ^ and  which  appear  in  the  momentm  and  energy  equa¬ 

tions.  These  modifying  terms  depend  upon  the  unsteady  part  of  the  solu¬ 
tion. 

Another  characteristic  of  these  equations  is  related  to  the 
way  in  which  the  quantity  ^varies.  For  example,  if  the  fluid  far  from 
the  cylinder  is  vibrating  in  an  arbitrary  direction,  Eq.  (6)  represents 
the  xinsteady  component  of  the  flow  far  from  the  cylinder.  Thus,  utili*- 
ing  Lin's  notation: 

=  2  acj  sin  <^t  sin  (•-  +  X)  (24) 


=  2  —  sin  CJ  t  cos  (|  +  X) 

Substituting  Eqs.  (24)  and  (25)  into  the  defining 
time  averages  gives: 


(25) 

Eq.  (15)  and  taking  the 


=  al„  (|  *  X) 


cos 


(- 


X)  - 


(26) 
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V 


Inspection  of  Eq.  (26)  shows  that: 

At  y  =  0,  Ui  =  Vj.  =  0,  ?=  2  sin  (|  +X)  cos  (|  +>.)  (27) 

At  0,  ^^0 

The  preceding  calculation  points  out  the  interesting  fact  that  the 
quantity  J  has  its  maximum  effect  at  the  wall.  Also,  this  calculation 
shows  that  the  quantity  acJ  is  an  important  parameter  in  determining 
the  amoiant  by  which  the  flow  is  modified  due  to  the  oscillations. 

This  conclusion  agrees  with  experimental  results. 

Another  Important  characteristic  of  the  flow  which  Lin  demon¬ 
strated  by  transforming  Eq.  (15)  is  that  a  pressure  gradient  outside 
the  boundary  layer  is  necessary  in  order  that  the  quantity^ be  different 
from  zero. 

Lin's  method  is  an  iterative  method  of  solution  and  makes  use 
of  the  basic  equations  in  a  dimensional  form.  For  most  other  analytical 
methods  of  attack,  the  differential  equations  in  dimensionless  form  are 
the  most  useful.  These  dimensionless  equations  also  supply  important 
hints  as  to  which  dimensionless  groups  best  describe  the  problem.  Such 
dimensionless  equations  will  be  discussed  in  the  next  chapter. 


I 
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III.  DIMENSIONLESS  EQUATIONS  AND  SPECIFIC  SOLUTIONS 
3.1  Stream  Function 


A  simplification  of  the  differential  equations  and  boundary 
conditions  can  be  obtained  by  use  of  the  stream  function,  ^ ,  The 
continuity  equation  is  satisfied  identically  by  ;  thus, 

u  =  V  =  ~ 

Dy  '  ” 

The  energy  and  momentm  equations  take  the  following  form; 


(28) 


•s — ^  - r — ^  -  ■:? - r — =  2  a  w  cos  w  t  sin  (r  +  X  )  + 

dydt  9y  «)yc>x  3x  Q  2  '•R  ' 


4  -  sin*"  t  sin  ^  X  )  cos  (f  X  )  +0  +  g8  O  sin^  (29) 

»y 


at  dy  dx  ■  ^x  3y  ^Cp  ^y2 

B.C.  y  =  0,  .^  =  .^  =  0,  6=  e 

’  d  X  ^y  w 

oO  ,  2  a  sin  wt  sin  (j  X  ) ,  0  0 


(30) 

(31) 


In  deriving  the  above  expression  for  the  momentum  equation,  Eq.  (25),  the 
expression  for  the  free-stream  velocity,  Eq.  (6) ,  was  substituted  into 
the  non-steady  Bernonlli  equation,  Eq.  (4),  and  the  result  of  this 
substitution  was  used  to  eliminate  the  pressure  term  from  Eq.  (2) . 

3.2  General  Dimensionless  Equations 

Eqs.  (29)  and  (30)  and  the  boundary  conditions  (31)  can  be 
made  dimensionless  by  the  following  substitutions: 
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t  =  Cjt 


X  = 


M  = 


<!> 


a. 


(32) 


y  =  •Cy 

In  terms  of  these  dimensionless  quantities!  Eqs.  (29) »  (30) ,  and  (3l) 
take  the  follovdng  form; 

2  ^ 

+  Ma  Maa  -  Ma  Maa  =  -  cos  t  sin  (x  +  X)  + 

•(  ir  ^  y  yx  X  yy  (^y)2 


+  4  (^y )  sin  t  sin  (x  +  X)  cos  (x  +  X)  +  ^ 


Maa  A  + 

yyy 


(33) 


(o<r)^  ^ 


A 

sin  X 


A  A 

yy 


B.C.  y  =  0,  M*  =  M*  =  0,  ^  =  1 

y->  "O,  Ma  sin  t  sin  (x  +  X) ,  ^  "*0 


(34) 


(35) 


where  the  subscripts  denote  partial  differentiation  with  respect  to 
dimensionless  coordinates. 

3.3  Specific  Dimensionless  Equations 

From  this  point  on,  all  coordinates  in  this  report  will  be 
dimensionless,  and  no  flexes  will  be  used  to  designate  them  as  such, 
in  order  to  simplify  the  writing. 
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The  quantities  at  and  TT  »Aich  appear  in  Eqs.  (33)  through 
(35)  are  as  yet  unspecified.  These  quantities  cam  now  be  chosen  so 
as  to  make  as  mamy  of  the  coefficients  in  Eq.  (33)  equal  to  unity  as 
possible.  Two  is  the  maximum  number  that  cam  be  made  equal  to  unity 
slmultameously.  Four  possible  choices  for  «(  amd  iT  indicated  in 
Table  1. 


When  the  pairs  of  values  for  and  Jf  in  Table  1  are  substi¬ 
tuted  into  Eqs.  (33) *  (34) »  amd  (35) ,  the  following  four  sets  of  equa¬ 
tions  are  obtained: 

Case  I, 

f<*»^  (1)  (1)  (1)  (1)  (1)  _  ,av  V  +  4.  4. 

^gOe  ^  )  cos  t  sin  (x  +  A.)  + 


yt  y  yx  X  yy 


2  7 

+  4  (|)  (^^Q  )  sin^  t  sin  (x  +  X)  cos  (x  +  X)  +  +  ^sin  x  (36) 


^g|3© 

^  t  *  ^  X  ■  "«' V  y  =  fc  ^ 

B.C.  y  =  0,  =0,  ^  =  1 


yyy 


yy 


y-»eo,  M^^^->2(^)  ("^-)  sin  t  sin  (x  +A),  -*0 


(37) 


(38) 


Case  II, 

(2)  „(2)  „(2)  „(2)  „(2)  _  „  a  +  .  /  ^  a  n  ^ 

V  ”y  V  “  X  V  R  t  sin  (x  +  A)  + 

2 

+  4  (f)  sln^  t  sin  (x  +  X  )  cos  (x  +  X)  +  + 

R  yyy 

*  ^  "1"  > 


(39) 
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(40) 


B.C.  y  =  0,  =0,  ^  =  1 

y->oe,  sin  t  sin  (x  +  X) , 

Case  III, 

a  yt  y  yx  X  "yy  ~  ^  ^  (*  +  >.)  + 

2 

+  4  sin  t  sin  (x  +  X)  cos  (x  +  X)  +  + 

yyy 

+  fi^^w  !>■}  A 

^  ^2  )  ^  sin  X 

'!>  =  <fi 

B.C.  y  =  0,  =0,  ^  =  1 

y-*0O,  2  sin  t  sin  (x  +  X),  ^-►O 

Case  IV, 

*J^  -  2  +  4  /• 

V  a  yt  y  yx  x  ^yy  '  2  cos  t  sin  (x  +  X)  + 

+  4  (p)  sin  t  sin  (x  +  X)  cos  (x  +  X)  +  + 

vtr»r 


yy 


g/9e/ 


,3 


yyy 


isin  X 


»>^  "'VSr’ur 


(41) 


(42) 


(43) 


(44) 


(45) 


yy 


(46) 
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(47) 


B.C.  y  =  0,  =0,  ^  =  1 

y-+  *0  ,  — *  2^^  sin  t  sin  (x  +  A  )  ,  ^  0 

The  nvunbers  in  parentheses  placed  above  the  dimensionless  stream 
function,  M,  in  the  preceding  sets  of  equations,  refers  to  the  Case  No. 
in  Table  1. 

3.4  Dimensionless  Parameters 

The  preceding  sets  of  equations  show  that  there  are  four  dimen¬ 
sionless  parameters  that  fully  describe  the  flow;  namely,  a  fluid  para¬ 
meter  represented  by  the  Prandtl  number,  a  temperature  parameter  represented 
by  the  Grashoff  number,  a  vibration  parameter  (vibration  parameter  A) 
represented  by  the  ratio  of  vibration  amplitude  to  cylinder  radius,  and 
another  vibration  parameter  (vibration  parameter  B)  which  may  be  repre¬ 
sented  by  any  of  four  different  dimensionless  groups.  The  four  sets  of 
dimensionless  parameters  corresponding  to  the  four  cases  in  Table  1  are 
listed  in  Table  2. 

3.5  Solutions  for  Uncoupled  Vibrations  and  Free  Convection 

Before  discussing  the  properties  of  each  of  the  four  sets  of 
equations  given  in  Section  3.3,  it  will  be  helpful  to  mention  here  solu¬ 
tions  obtained  for  boundary-layer  flow  for  the  case  of  vibrations  in  the 
absence  of  heat  transfer  and  for  free  convection  in  the  absence  of  vibra¬ 
tions.  The  first  of  these  cases  is  that  of  a  cylinder  oscillating  in  a 
fluid  otherwise  at  rest  with  no  temperature  difference  between  the  cylin¬ 
der  and  the  fluid.  The  solution  to  this  problem  was  obtained  by 
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Schlichtlng^^^^  by  the  uee  of  an  Iterative  method.  Schllchting* s  choice 
of  a  dimensionlese  y-Ksoordlnate  was  based  on  an  "AC"  boundary-layer  thick¬ 
ness  as  indicated  in  Case  II. 

The  problem  of  natural  convection  without  vibrations  was  solved 
(23^ 

analytically  by  Hermann'  '  for  a  heated  horizontal  cylinder.  Hermann 
nondlmensionallzed  his  equations  by  the  same  substitutions  as  indicated 
in  Case  I.  The  equations  were  then  transformed  into  ordinary  differential 
equations  by  an  approximate  similarity  variable  and  solved  numerically. 

3,6  Complete  Equations i  Case  I 

Using  the  substitution  of  variables  for  Case  I,  Table  II, 

(2/) 

T.  Chlang'  ,  following  some  ideas  suggested  by  R.  J.  Nickerson,  attacked 
the  complete  analytical  problem  for  simultaneous  heat  transfer  and  verti¬ 
cal  vibrations  ( X  =  0) .  By  a  perturbation  method  of  solution,  Chlang 
obtained  numerical  values  that  satisfied  his  zeroth-order  equations;  he 
also  worked  out  the  form  for  the  higher-orde-^  equations.  The  numerical 
values  for  the  zeroth-order  equations,  which  are  the  same  equations  as 
for  steady  natural  convection,  agree  very  well  with  the  values  obtained 
by  Hermann.  T.  Chiang  used  a  different’ similarity  variable  than  Hemann, 
in  order  to  be  able  to  use  the  same  variable  for  the  higher-order  equa¬ 
tions.  A  short  siimmary  of  Chiang* s  perturbation  method  and  results  fol¬ 
lows.  For  more  complete  details,  the  reader  is  referred  to  the  original 
work^^^^ . 


In  Chiang's  approach,  the  dimensionless  stream  function, 
and  temperature  function,  ^ ,  are  expanded  in  power  series  in  terms  of  a./R: 


M 


(x,  y) 


(|)  (x,  y. 


t)  + 


(-) 


M, 


(1) 


(x,  y,  t)  + 


(48) 
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2 

^  (xi  y)  +  (x»  yt  t)  +  (|)  ^2^2  y*  ■*■ 


•  •  •  • 


(49) 


These  expressions  are  substituted  into  Eqs.  (36) «  (37) ,  and  (38) •  Then, 
by  equating  coefficients  of  like  powers  of  the  perturbation  parameter, 
a/ll,  the  following  set  of  equations  and  boundary  conditions  are  obtained: 
Zeroth-Order  Equations 


^  sin 

oy  oyx  ox  oyy  oyyy  •  o 


(L  -  (h  =  4-  <2^ 

oy  ox  or.  r  oy  Pr  n 


oyy 


B.C.  y  =  0,  =  0,  C/>  =  1 

'  'ox  oy  '  "  o 


►  ,  M 


(1) 


oy 


0,  ^  -0 


First-Order  Equations 


(50) 

(51) 


(52) 


oy  lyx 

ly  oyx  ox  lyy 

-  = 

lx  oyy 

(53) 

2  cos 

t  sin  X  + 

+  (£ 

oy  IT 

.  +  (A  -  (A 

lx  ly  or.  or  r 

ly  '  ^x^  ^oj  ~  Pr  J 

(54) 

B.C.  y  = 

0,  M<«  = 

:K<«  =0,  ^j=0 

(55) 

y-e 

-♦  2  ^  sin  t  sin  x,  ^  ^ 

-*•  0 

Secord-Order  Equations 

^  ^yt  oy  2yx  ly  lyx  2y  oyx  ox  2yy  T.x  T.yy 
“  ”2x^  ^oyy  ^  ^  ^  x  cos  x  +  sin  x 


(56) 
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^2y  '  «1*  ^ly  * 

-'4i'(^oy=fe  '^2yy 

B.C.  y  =  0,  «^>  =  =  0,  (^2  =  0 

y-**o  ,  0,  ^2“^  ° 

A  reduction  of  the  zeroth-order  equations  to  ordinary  differen¬ 
tial  equations  can  be  achieved  by  the  follovdng  substitution: 

=  X  (y)  +  (y)  +  F2  (y)  (59) 

0  o  =  (y)  +  X^  (y)  +  X^  G2  (y)  (60) 


The  values  for  the  functions  F^,  F^,  F^*  G^,  G^^,  G2,  and  their  derivatives 
have  been  calculated  and  are  given  in  Tables  3»  4>  and  5  of  the  Appendix. 

From  Eqs.  (32),  (59),  (60)  and  Tables  3  to  5,  it  can  be  shown 
that  the  local  Nusselt  number  and  the  dimensionless  velocity  and  tempera¬ 
ture  are,  respectively, 


Nu 


0.3702  +  0.0161 


0.00007 


(61) 


X  F  +  X^  F,  +  X^  F„ 
oy  ly  2y 

=  G^  +  X^  G^  +  X^  G2 


(62) 


(63) 
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By  using  exptuisions  similar  to  Eqs.  (59)  and  (60),  the  higher- 
order  equations  can  also  be  reduced  to  ordinary  differential  equations. 

An  advantage  of  this  perturbation  method  is  that  all  dimen¬ 
sionless  parameters,  except  the  Prandtl  number,  have  been  eliminated 
from  the  zeroth-order  equations,  and  hence  these  equations  need  to  be 
solved  only  once  for  each  Prandtl  number.  The  higher-order  equations 
are  linear;  this  fact  constitutes  another  advantage  of  the  perturbation 
method. 

Two  important  disadvantages  are  present  in  this  method.  The 
first  concerns  the  question  of  how  to  relate  the  theoretical  solution 
to  the  experimental  results.  Experiments  indicate  that  the  quantity  &<•> 
is  an  important  factor  for  determining  the  nature  of  the  flow.  However, 
this  combination  does  not  appear  anywhere  among  the  dimensionless  parame¬ 
ters  of  this  case  (Case  I,  Table  2).  This  difficulty  might  prove  imagi¬ 
nary  when  the  equations  are  solved  and  the  numerical  results  are  compared 
with  the  experimental  results.  However,  the  second  disadvantage  is 
very  real  and  concerns  the  dimensionless  parameter,  ,  which  appears 
in  the  higher-order  equations  (in  addition  to  the  Prandtl  number).  Since 
this  frequency  parameter,  a,  contains  the  cylinder  temperature  and  the 
frequency  of  vibration,  it  necessitates  a  completely  new  solution  of  the 
higher-order  equations  for  almost  every  different  value  of  cylinder 
temperature  and  frequency. 

These  disadvantaiges  provided  the  motivation  to  investigate 
the  possibility  of  using  other  sets  of  parameters  than  those  in  Case  I. 


These  other  sets  of  parameters-— Cases  II,  III,  and  IV  in  Table  2--over- 
come  some  of  the  above  disadvantages,  but  new  disadvantages  appear  In 
their  place.  However,  by  studying  these  additional  cases,  a  greater 
insight  into  the  problem  is  obtained. 

3.7  Complete  Equations,  Case  II 

This  set  of  equations  Is  very  similar  to  those  of  Case  I  and 


can  be  expanded  in  a  power  series  like  Eqs.  (48)  and  (49) .  When  this 
Is  done,  the  following  sets  of  equations  are  obtained: 

Zeroth-Order  Equations 


- 

= 

+  Gr  (i^)  1 

®  sin  X 

(64) 

oy 

oyx 

ox 

oyy 

oyyy 

r  0 

M^2) 

6  - 

s  = 

h  (t> 

1 

(65) 

oy 

T  ox 

ox 

r  oy 

Pr  T 

oyy 

B.C. 

o 

It 

>. 

oy 

II 

=  0, 

=  1 

y->oo. 

oy 

o 

t 

-♦0 

(66) 

First-Order  Equations 

2 

+  +  Gr  (|i)  ^  ^  sin  X  (67) 

It  -  01X  ^  ^y  4>o.  -  <t>^,  -  <t>o,  Tr  <Plyy  ^^S) 

B.C.  y  =  0,  =  Mjj)  =0,  4>^  =  0 

y  00  ,  -*  2  sin  t  sin  (x  +  ^) ,  ^  ^  0 


(69) 
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Second-Order  Equations 

(2)  .  „(2)  „(2)  ,  „(2)  ^(2)  .  J2)  „(2)  „(2)  „(2)  _  w(2)  J2) 

^yt  ^  ”oy  ”2yx  ^  ^y  ”lyx  ^  ”2y  ”oyx  "  ”ox  ^  ”lx  ^yy  ^ 

=  4  sln^  t  cos  (x  +  X)  sin  (x  +  X)  +  + 


'Am  a  a 

2x  oyy 
2 


iyyy 


(70) 


sin  X 


+  Gr(|)  <!> 

Ac  «o?  <Z^2x  "  <^lx  '  >4^’  0..  -  ««  02y  -  A, 

-"Sx’  0oy  =  P;  02yy 
B.C.  y  =  0,  =0,  02"  ° 

y-*  oo,  0,  02'^^ 

An  interesting  possibility  presented  by  this  set  of  equations 

is  to  try  to  eliminate  the  dimensionless  parameters  from  the  first  and 

second-order  equations.  This  could  be  done  by  absorbing  the  quantity 

Gr  into  a  new  definition  of  the  temperature  function.  Since  the 

energy  equation  is  linear  and  the  temperature  boundary  conditions  are 

homogeneous,  this  step  presents  no  great  difficulty. 

However,  one  difficulty  still  remains  with  this  set  of  equa- 

t  2 

tions:  Since  the  quantity  Gr  (~)  appears  in  the  momentum  equation, 

Eq.  (64),  the  zeroth-order  equations  have  to  be  solved  for  each  new 
condition  of  temperature  or  frequency.  Eliminating  this  quantity  from 
Eq.  (64),  by  defining  a  new  temperature  function,  does  not  overcome  this 


(71) 


(72) 
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difficulty  because  the  temperature  boundary  conditions,  Eq.  (66),  are 

not  homogeneous.  Hence,  although  the  differential  equation  would 

remain  unchanged,  the  boxmdary  conditions  at  the  cylinder  wall  would 
j  r  2 

become  (p ^  -  Gr  (^)  instead  of  being  equal  to  unity.  This  boundary 
condition  would  change  for  each  new  wall  temperature  or  frequency, 
thereby  necessitating  a  complete  numerical  solution  for  each  new  condi¬ 
tion. 

There  is  a  possible  way  of  circumventing  the  main  disadvan¬ 
tages  of  Cases  I  and  II.  This  would  be  to  combine  the  zeroth-order 
equations  of  Case  I  with  the  higher-order  equations  of  Case  II.  Then, 
the  equations  would  have  to  be  solved  only  once  for  each  Prandtl  number. 
This  possibility  should  be  investigated  in  the  future. 

The  advantage  of  linear  higher-order  equations  is  present  in 
Case  II,  as  in  Case  I,  Also,  as  in  Case  I,  the  quantity  acj  does  not 
appear  anywhere  among  the  dimensionless  parameters  of  Case  II — for  rea¬ 
sons  stated  in  section  3.6,  the  nonappearance  of  the  quantity  aU  may  be 
a  disadvantage  in  solving  the  problem. 

3.8  Complete  Equations,  Case  III 

For  Case  III,  expanding  the  equations  by  perturbation  methods 
does  not  give  physically  meaningful  equations,  since  the  highest  deriva¬ 
tive  terms  do  not  appear  in  the  zeroth-order  solution.  Therefore,  other 

mathematical  approaches  must  be  tried.  Integral  approaches  similar  to 

(19) 

those  used  by  Llghthill  might  prove  successful. 

In  this  case,  unlike  the  two  previous  cases,  the  product  a<j 
does  appear  in  the  dimensionless  parameters.  It  appears  in  vibration 
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parameter  B,  the  vibrational  Reynolds  number.  Experiments  indicate  that 
this  product,  a<*>  ,  is  an  important  characteristic  of  the  flow.  For  this 
reason,  this  form  of  the  equations  might  bring  fruitful  results  with  addi¬ 
tional  investigation. 

3.9  Complete  Equations,  Case  IV 

Case  IV  appears  to  have  no  special  advantages  over  Case  III 
which  it  resembles.  The  appearance  of  fractional  powers  of  the  ratio 
a/R  is  an  additional  complication  without  compensating  advantages. 


25 


IV.  SUMMARY 


In  this  report,  the  basic  equations  and  the  boundary  conditions 
which  govern  the  problem  of  coupled  transverse  vibrations  and  free  convec¬ 
tion  from  a  heated  horizontal  cylinder  are  presented.  By  applying  a 

(21^ 

method  developed  by  C.  C.  Lin'  '  some  important  characteristics  of  these 
equations  are  demonstrated.  It  la  shown  that  the  steady-flow  solution  is 
modified  by  the  presence  of  harmonic  oscillations  only  when  pressure  gradi¬ 
ents  occur.  These  modifications  have  their  most  pronotinced  effect  on  the 
fluid  closest  to  the  surface.  Furthermore,  Lin's  method  indicates  that 
the  product  aCJ  is  an  important  parameter  for  determining  the  magnitude 
of  the  changes  produced  in  the  steady-flow  solution  by  the  presence  of 
harmonic  oscillations;  this  fact  agrees  with  experimental  observations. 

Lin's  method  is  an  iterative  method  of  solution  and  makes  use 
of  the  basic  equations  in  a  dimensional  form.  For  most  other  analytical 
methods  of  attack,  the  differential  equations  in  dimensionless  form  are 
the  most  useful.  For  this  reason,  the  basic  differential  equations  were 
converted  into  dimensionless  fom.  An  examination  of  these  dimensionless 
equations  shows  that  four  dimensionless  parameters  are  needed  to  fully 
describe  the  flow.  Table  2  gives  four  possible  sets  of  dimensionless 
parameters;  the  differential  equations  associated  with  these  four  sets 
of  parameters  are  presented  in  Chapter  III. 

Considerable  work  has  been  done  on  the  first  set  of  parameters 
in  Table  2  (Case  l)  by  T.  Chiang  .  The  main  advantage  of  formulating 
the  problem  in  terms  of  the  parameters  in  Case  I  is  the  simplicity  of 
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the  zeroth-order  equations  obtained  when  a  perturbation  expansion  in 
terms  of  a/R  is  applied.  These  zeroth-order  eq\xations,  which  are  identi¬ 
cal  with  the  equations  for  free  convection  without  vibration,  have  been 

solved  by  Chiang,  and  the  results  of  this  solution  agree  very  well  with 

(23) 

a  solution  obtained  by  Hermann'  '  who  used  other  methods. 

For  the  dimensionless  parameters  in  Case  II,  the  higher-order 
equations,  obtained  by  applying  a  perturbation  expansion  in  terms  of  a/E., 
are  simpler  them  the  corresponding  higher-order  equations  for  Case  I. 

By  combining  the  perturbation  equations  for  Case  I  and  Case  II,  the 
advantages  of  each  might  be  retained  for  a  complete  solution.  This 
possibility  merits  further  investigation. 

The  dimensionless  parameters  for  Case  III  contain  the  product 
aco  ,  which  experiment  has  shown  to  be  a  controlling  variable;  for  this 
reason.  Case  III  represents  a  promising  line  of  attack.  Case  IV  seems 
to  present  no  particular  advantages  over  the  first  three  cases  in  obtain¬ 
ing  a  solution  to  the  problem. 
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TABLE  1 


Values  for  cH  and  Tf  used  to  nondimensionalize  the  basic  differential  equations 


TABLE  2 

Sets  of  dimensionless  parameters 


Case  No. 

I 

II 

III 

IV 

Fluid  Parameter 

Pr 

Pr 

Pr 

Pr 

Temperature  Parameter 

Gr 

Gr 

Gr 

Gr 

Vibration  Parameter  (a) 

a 

R 

a 

R 

a 

R 

a 

R 

Vibration  Parameter  (B) 

2  r  uhi 

g{3ew 

-2. 

(Rel  -  ^ 

^af  ^ 

Va?<jR 

)> 
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APPENDIX 
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TABLE  3 

Tabulation  of  the  fvinctions  F  and  G  and  their  derivatives 

o  o 

(Pr  =0.70) 


■■1 

F 

oyy 

F 

oy 

F 

0 

G 

oy 

G 

0 

0.0 

0.8595 

0.0000 

0.0000 

-0.3702 

1.0000 

.1 

.7614 

.0810 

.0041 

-.3702 

.9630 

.2 

.6677 

.1524 

.0159 

-.3699 

.9260 

.3 

.5789 

.2147 

.0343 

-.3693 

.8890 

.4 

.4953 

.2684 

.0585 

-.3681 

.8521 

.5 

.4171 

.3140 

.0877 

-.3662 

.8154 

.6 

.3446 

.3520 

.1211 

-.3636 

.7789 

.7 

.2778 

.3831 

.1579 

-.3600 

.7427 

.8 

.2166 

.4078 

.1975 

-.3556 

.7069 

.9 

.1611 

.4266 

.2393 

-.3502 

.6716 

1.0 

.1109 

.4401 

.2826 

-.3439 

.6369 

1.1 

.0661 

.4490 

.3271 

-.3366 

.6029 

1.2 

.0263 

.4535 

.3723 

-.3285 

.5696 

1.3 

-.0088 

.45U 

.4177 

-.3195 

.5372 

1.4 

-.0393 

.4519 

.4630 

-.3098 

.5057 

1.5 

-.0655 

.4467 

.5080 

-.2995 

.4753 

1.6 

-.0878 

.4390 

.5523 

-.2886 

.U59 

1.7 

-.1066 

.4292 

.5957 

-.2772 

.4176 

1.8 

-.1219 

.4178 

.6381 

-.2655 

.3904 

1.9 

-.1343 

.4049 

.6792 

-.2535 

.3645 

2.0 

-.U39 

.3910 

.7190 

-.2414 

.3397 

2.1 

-.1511 

.3762 

.7574 

-.2292 

.3162 

2.2 

-.1560 

.3609 

.7943 

-.2171 

.2939 

2.3 

-.1591 

.3451 

.8296 

-.2051 

.2728 

2.4 

-.1605 

.3291 

.8633 

-.1933 

.2529 

2.5 

_»1604 

.3130 

.8954 

-.1818 

.2341 

2.6 

-.1590 

.2971 

.9259 

-.1703 

.2165 

2.7 

-.1566 

.2813 

.9543 

-.1597 

.2000 

2.8 

-.1534 

.2658 

.9821 

-.U92 

.1846 

2.9 

-.1493 

.2506 

1.0080 

-.1392 

.1701 

3.0 

-.1448 

.2359 

1.0320 

-.1296 

.1567 

3.2 

-.1343 

.2080 

1.0770 

-.1118 

.1326 

3.4 

-.1229 

.1822 

1.1160 

-.0959 

.1119 

3.6 

-.1111 

.1588 

1.1500 

-.0818 

.0941 

3.8 

-.0995 

.1378 

1.1790 

-.0695 

.0790 

4.0 

-.0883 

.1190 

1.2050 

-.0588 

.0662 

-.0778 

-.0681 

.1024 

.0879 

1.2270 

1.2460 

-.0496 

-.0417 

.0554 

.0463 

-.0593 

.0751 

1.2620 

-.0350 

.0386 

4.8 

-.0514 

.0641 

1.2760 

-.0293 

.0322 

5.0 

-.0U3 

.0545 

1.2880 

-.0245 

.0269 
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TABLE  3  (cont'd) 


(Pr  =  0.70) 


y 

F 

oyy 

F 

oy 

F 

o 

G 

oy 

G 

0 

5.2 

-.0381 

.0463 

1.2980 

-.0204 

.0224 

5.4 

-.0326 

.0392 

1.3070 

-.0170 

.0186 

5.6 

-.0278 

.0332 

1.3140 

-.0142 

.0155 

5.8 

-.0237 

.0281 

1.3200 

-.0118 

.0129 

6.0 

-.0202 

.0237 

1.3250 

-.0098 

.0108 

6.5 

-.0133 

.0154 

1.3350 

-.0061 

.0069 

7.0 

-.0087 

.0100 

1.3410 

-.0038 

.0044 

7.5 

-.0057 

.0064 

1.3450 

-.0024 

.0029 

8.0 

-.0037 

.0041 

1.3480 

-.0015 

.0019 

8.5 

-.0025 

.0026 

1.3490 

-.0009 

.0013 

9.0 

-.0017 

.0016 

1. 3500 

-.0006 

.0009 

9.5 

-.0011 

.0009 

1.3510 

-.0004 

.0007 

10.0 

-.0008 

.0004 

1.3510 

-.0002 

.0006 

10.5 

-.0006 

.0000 

1.3510 

-.0001 

.0005 
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TABLE  4 


Tabulation  of  the  functions  and  and  their  derivatives 


(Pr  =0.70) 


■ 

^lyy 

ly 

"l 

ly 

°1 

0.0 

-0.09342 

0.0000 

0.0000 

0.01609 

0.0000 

.1 

-.0772 

-.0085 

-.0004 

.0161 

.0016 

.2 

-.0620 

-.0155 

-.0017 

.0160 

.0032 

.3 

-.0481 

-.0210 

-.0035 

.0159 

.0048 

.4 

-.0354 

-.0251 

-.0058 

.0157 

.0064 

.5 

-.0242 

-.0281 

-.0085 

.0154 

.0080 

.6 

-.0143 

-.0300 

-.0114 

.0149 

.0095 

.7 

-.0058 

-.0310 

-.0144 

.0143 

.0109 

.8 

.0013 

-.0312 

-.0176 

.0136 

.0123 

.9 

.0073 

-.0308 

-.0207 

.0127 

.0136 

1.0 

.0122 

-.0298 

-.0237 

.0118 

.0149 

1.1 

.0160 

-.0284 

-.0266 

.0107 

.0160 

1.2 

.0189 

-.0266 

-.0294 

.0095 

.0170 

1.3 

.0210 

-.0246 

-.0319 

.0083 

.0179 

1.4 

.0224 

-.0224 

-.0343 

.0070 

.0187 

1.5 

.0232 

-.0202 

-.0364 

.0057 

.0193 

1.6 

.0234 

-.0178 

-.0383 

.0045 

.0198 

1.7 

.0233 

-.0155 

-.0400 

.0032 

.0202 

1.8 

.0227 

-.0132 

-.0414 

.0020 

.0205 

1.9 

.0219 

-.0110 

-.0426 

.0008 

.0206 

2.0 

.0208 

-.0088 

-.0436 

-.0003 

.0206 

2.2 

.0182 

-.0049 

-.0450 

-.0023 

.0204 

2.4 

.0153 

-.0016 

-.0456 

-.0039 

.0197 

2.6 

.0124 

.0012 

-.0456 

-.0051 

.0188 

2.8 

.0096 

.0034 

-.0452 

-.0060 

.0177 

3.0 

.0070 

.0051 

-.0443 

-.0065 

.0165 

3.2 

.0048 

.0062 

-.0432 

-.0068 

.0151 

3.4 

.0029 

.0070 

-.0418 

-.0068 

.0138 

3.6 

.0013 

.0074 

-.0404 

-.0066 

.0124 

3.8 

.0001 

.0075 

-.0389 

-.0063 

.0111 

-.0009 

.0074 

-.0374 

-.0060 

.0099 

-.0016 

.0072 

-.0359 

-.0055 

.0088 

-.0021 

.0068 

-.0345 

-.0051 

.0077 

4.6 

-.0024 

.0064 

-.0332 

-.0046 

.0067 

4.8 

-.0026 

.0058 

-.0320 

-.0041 

.0059 

5.0 

-.0027 

.0053 

-.0308 

-.0037 

.0051 

5.5 

-.0025 

.0040 

-.0286 

-.0027 

.0035 

6.0 

-.0021 

.0028 

-.0269 

-.0020 

.0023 

6.5 

-.0016 

.0019 

-.0257 

-.0014 

.0015 

7.0 

-.0012 

.0012 

-.0249 

-.0010 

.0009 

7.5 

-.0008 

.0007 

-.0244 

-.0007 

.0005 

8.0 

-.0005 

.0004 

-.0242 

-.0004 

.0002 

8.5 

-.0003 

.0002 

-.0241 

-.0003 

.0001 

8.7 

-.0002 

.0001 

-.0240 

-.0002 

.0000 

TABLE  5 


Tabulation  of  the  functions  F2  and  G2  and  their  derivatives 


(Pr  =0.70) 


y 

^2yy 

2y 

^2 

°2 

0.0 

0.00245 

0.00000 

0.00000 

-0.00007 

0.00000 

.1 

.00165 

.00001 

-.00007 

-.00001 

.2 

.00094 

.00004 

-.00007 

-.00001 

.3 

.00034 

.00008 

-.00007 

-.00002 

.4 

-.00015 

.00012 

-.00006 

-.00003 

.5 

-.00053 

.00016 

-.00005 

-.00003 

•6 

-.00081 

.00019 

-.00003 

-.00004 

.7 

-.00099 

.00021 

.00022 

-.00001 

-.00004 

-.00108 

.00011 

.00023 

.00001 

-.00004 

-.00111 

-.00000 

.00024 

.00003 

-.00004 

1.0 

-.00107 

-.00011 

.00023 

.00005 

-.00003 

-.00099 

-.00021 

.00021 

.00007 

-.00003 

-.00087 

-.00031 

.00019 

.00008 

-.00002 

-.00072 

-.00039 

.00015 

.00009 

-.00001 

-.00055 

-.00045 

.00011 

.00009 

-.00000 

-.00036 

-.00050 

.00006 

.00009 

.00001 

-.00017 

-.00052 

.00001 

.00008 

.00001 

.00002 

-.00053 

-.00004 

.00006 

.00002 

1  1.8 

.00022 

-.00052 

-.00009 

.00004 

.00003 

MS 

,00042 

-.00049 

-.00014 

.00001 

.00003 

hM 

.00061 

-.00044 

-.00019 

-.00003 

.00003 

.00080 

-.00036 

-.00023 

-.00007 

.00002 

.00099 

-.00027 

-.00026 

-.00012 

.00001 

PB9 

.00118 

-.00017 

-.00028 

-.00017 

.00000 
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